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Abstract 



We discuss Lambda polarization in semi-inclusive proton-proton collisions, 
with one of the protons longitudinally polarized. The hyperfine interaction 
responsible for the A-iV and E-A mass splittings gives rise to flavor asym- 
metric fragmentation functions and to sizable polarized non-strange fragmen- 
tation functions. We predict large positive Lambda polarization in polarized 
proton-proton collisions at large rapidities of the produced Lambda, while 
other models, based on SU(3) flavor symmetric fragmentation functions, pre- 
dict zero or negative Lambda polarization. The effect of S° and E* decays is 
also discussed. Forthcoming experiments at RHIC will be able to differentiate 
between these predictions. 
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I. INTRODUCTION 



Measurements of the polarization dependent structure function, gi, in deep inelastic 
scattering [|TJ have inspired considerable experimental and theoretical effort to understand 
the spin structure of baryons. While most of these studies concern the spin structure of 
the nucleons, it has become clear that similar measurements involving other baryons would 
provide helpful, complementary information The Lambda baryon plays a special role 

in this respect. It is an ideal testing ground for spin studies since it has a rather simple 
spin structure in the naive quark parton model. Furthermore, its self-analyzing decay makes 
polarization measurements experimentally feasible. 

Forthcoming experiments at RHIC could measure the polarization of Lambda hyperons 
produced in proton-proton collisions with one of the protons longitudinally polarized, p^p — ► 
A^X. The polarization dependent fragmentation function of quarks and gluons into Lambda 
hyperons can be extracted from such experiments. These fragmentation functions contain 
information on how the spin of the polarized quarks and gluons is transferred to the final 
state Lambda. The advantage of proton proton collisions, as opposed to e + e~ annihilation, 
where A production and polarization is dominated by strange quark fragmentation, is that 
Lambdas at large positive rapidity are mainly fragmentation products of up and down 
valence quarks of the polarized projectile. Thus, the important question, intimately related 
to our understanding of the spin structure of baryons, of whether polarized up and down 
quarks can transfer polarization to the Lambda can be tested at RHIC ||. 

In a previous publication, we have shown that the hyperfine interaction, responsible for 
the A-N and S°-A mass splittings leads to non-zero polarized non-strange quark fragmen- 
tation functions [[TJJ]. These non-zero polarized up and down quark fragmentation functions 
give rise to sizeable positive A polarization in experiments where the strange quark fragmen- 
tation is suppressed. On the other hand, predictions based either on the naive quark model 
or on SU(3) flavor symmetry predict zero or negative Lambda polarization 0. 

In section II, we briefly discuss fragmentation functions and show how the hyperfine 
interaction leads to polarized non-strange fragmentation functions. We fix the parameters 
of the model by fitting the data on A production in e + e~ annihilation. In section III, we 
discuss A production in pp collisons at RHIC energies. We point out that the production 
of A's at high rapidities is dominated by the fragmentation of valence up and down quarks 
of the polarized projectile, and is ideally suited to test whether non-strange quarks transfer 
their polarization to the final state A. We predict significant positive A-polarization at large 
rapidities of the produced A. 



II. FRAGMENTATION FUNCTIONS 



Fragmentation functions can be defined as light-cone Fourier transforms of matrix ele- 
ments of quark operators [ |T^ . |T3| 

\ D U*) = \Y,j f r e- jP+r/z Tr{r(0|^(0)|A(P5);n( Pn ))(A(P5);n( Pn )|^(r)|0)}, (1) 

where, V is the appropriate Dirac matrix; P and p n refer to the momentum of the produced 
A and of the intermediate system n; S is the spin of the Lambda and the plus projections 
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of the momenta are defined by P + = (P° + P 3 ). z is the plus momentum fraction of the 
quark carried by the produced A. 

Translating the matrix elements, using the integral representation of the delta function 
and projecting out the light-cone plus and helicity ± components we obtain 



^ A (*) = ^]p[(i/*-i)p + -^ (2) 



Here, ip+ = ~7_7 + ~(l ± 75)^, and we have defined j± = 4g(7o ± 73). The fragmentation 
function of an anti-quark into a A is given by Eq. (|j), with ip + replaced by 

-Df A (z) = ^Y,m/z-±)P + -PmO\^ ± mMPS ll );n(p n ))\\ (3) 

Z Z\fZ n 



D^ A can be interpreted as the probability that a quark with positive/negative helicity frag- 
ments into a A with positive helicity and similarly for antiquarks. 

The operator ip + (^+) either destroys a quark (an antiquark) or it creates an antiquark 
(quark) when acting on the A on the right hand side in the matrix elements. Thus, whereas, 
in the case of quark fragmentation, the intermediate state can be either an anti-diquark 
state, qq, or a four- quark- antiquark state, qqqq, in the case of antiquark fragmentation, only 
four- antiquark states, qqqq, are possible assuming that there are no antiquarks in the A. 
(Production of A's through coupling to higher Fock states of the A is more complicated and 
involves higher number of quarks in the intermediate states. As a result it would lead to 
contributions at lower z values.) Thus, we have 

• (la) q — > qqq + qq = A + qq 

• (lb) q — > qqq + qqqq = A + qqqq, 
for the quark fragmentation and 

• (2) q — > qqq + qqqq = A + qqqq, 

for the antiquark fragmentation. 

While, in case (la), the initial fragmenting quark is contained in the produced Lambda, 
in case (lb) and (2), the Lambda is mainly produced by quarks created in the fragmentation 
process. Therefore, we not only expect that Lambdas produced through (la) usually have 
larger momenta than those produced through (lb) or (2) but also that Lambdas produced 
through (la) are much more sensitive to the flavor spin quantum numbers of the fragment- 
ing quark than those produced through (lb) and (2). In the following we assume that (lb) 
and (2) lead to approximately the same fragmentation functions. In this case, the differ- 
ence, D q A — Dq\, responsible for leading particle production, is given by the fragmentation 
functions associated with process (la). 

Similar observations also follow from energy-momentum conservation built in Eqs. (0) 
and (0). The delta function implies that the function, D q (z)/z, peaks at |14 



(4) 

M + M n K J 
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Here, M and M n are the mass of the produced particle and the produced system, n, and we 
work in the rest frame of the produced particle. We see that the location of the maxima of 
the fragmentation function depends on the mass of the system n. While the high z region is 
dominated by the fragmentation of a quark into the final particle and a small mass system, 
large mass systems contribute to the fragmentation at lower z values. The maxima of the 
fragmentation functions from process (la) are given by the mass of the intermediate diquark 
state and that of the the fragmentation functions from the processes (lb) and (2b) by the 
masses of intermediate four quark states. Thus, the contribution from process (la) is harder 
than those from (lb) and (2). 

Energy-momentum conservation also requires that the fragmentation functions are not 
flavor symmetric. While the assertion of isospin symmetry, D u \ = D^a, is well justified, 
SU(3) flavor symmetry is broken not only by the strange quark mass but also by the hyperfine 
interaction. Let us discuss the fragmentation of a u (or d) quark and that of an s quark into 
a Lambda through process (la). While the intermediate diquark state is always a scalar 
in the strange quark fragmentation, it can be either a vector or a scalar diquark in the 
fragmentation of the non-strange quarks. The masses of the scalar and vector non-strange 



diquarks follow from the mass difference between the nucleon and the Delta ||10|| , while those 
of the scalar and vector diquark containing a strange quark can be deduced from the mass 
difference between £ and A ]ll]|. They are roughly m s ~ 650 MeV and m v ~ 850 MeV 



for the scalar and vector non-strange diquarks, and m' s « 890 MeV and m' v ~ 1010 MeV 
for scalar and vector diquarks containing strange quarks, respectively []1T,14|. According to 



Eq. (f|), these numbers lead to soft up and down quark fragmentation functions and to hard 
strange quark fragmentation functions. 

Energy-momentum conservation, together with the splitting of vector and scalar diquark 
masses, has the further important consequence that polarized non-strange quarks can trans- 
fer polarization to the final state Lambda. To see this we note that the probabilities for the 
intermediate state to be a scalar or vector diquark state in the fragmentation of an up or 
down quark with parallel or anti-parallel spin to the spin of the Lambda can be obtained 
from the SU(6) wave function of the A 

A T = — =[2s T Mo, + V2d^(us)i t i — d\us)i t0 + d^(us) 00 + 
2y 3 

-V2u l (ds)i )1 + u l (ds)i )0 -u l (ds) 0) o} . (5) 

While the u or d quarks with spin anti-parallel to the spin of the A are always associated 
with a vector diquark, u and d quarks with parallel spin have equal probabilities to be 
accompanied by a vector or scalar diquark. The fragmentation functions of non-strange 
quarks with spin parallel to the A spin are harder than the corresponding fragmentation 
functions with anti-parallel spins. Thus, AD u a is positive for large z values and negative 
for small z. Their total contribution to polarized Lambda production might be zero or very 
small. Nevertheless, AD U \ and AD^a can be sizable for large z values, since both D u a and 
AD U \ are dominated by the spin-zero component in the large z limit. Furthermore, they 
will dominate polarized Lambda production whenever the production from strange quarks 
is suppressed. 

The matrix elements can be calculated using model wave functions at the scale relevant to 
the specific model and the resulting fragmentation functions can be evolved to a higher scale 
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to compare them to experiments. In a previous paper , we calculated the fragmentation 
functions in the MIT bag model and showed that the resulting fragmentation functions 
give a very reasonable description of the data in e + e~ annihilation. Since the mass of the 
intermediate states containing more than two quarks are not known we only calculate the 
contributions of the diquark intermediate states in the bag model. The other contributions 
have been determined by performing a global fit to the e + e~ data. For this, we used the 
simple functional form 

D qA (z) =N q z a {l-zY (6) 

to parameterize D qA = D nA = D^ A = ...Di A and also set D gA = at the initial scale, 
H = 0.25 GeV. 

The fragmentation functions have to be evolved to the scale of the experiment, /i. The 
evolution of the non-singlet fragmentation functions in LO is given by |15| , |16l 



d /"l d? 7 

[D qA - D qA ](z^ 2 ) = / -P qq {-)[D qA - D qA ){z\y?). (7) 

J z Z Z 



d\n.jj? 

The singlet evolution equations are 

1 ±- 2 Y^D qA {z^ 2 ) = J^[P qq ^)^D qA (z,^) + 2n f P 9 ^)D gA (z'^ 2 )} 

rj rl A ] 2 7 7 

— ZV(^ 2 ) = / — [P^llDAz^ 2 ) + P gg (-)D gA (z',^)], (8) 

where the splitting functions are the same as those for the evolution of quark distributions. 
rif is the number of flavors. We used the evolution package of Ref. [O suitably modified 
for the evolution of fragmentation functions (interchanging the off-diagonal elements in the 
singlet case). 

The results of the LO fit to e + e~ data |18|-Q are shown in Fig. 1. The parameters of our 



fits are given in Table I. The bag model calculations for D qA — D qA and AD sA — AD sA were 
parametrized using the functional form of Eq. The fragmentation functions, AD uA — 
AD uA = AD dA — ADj A , change sign at some value of z, hence we parametrized them using 
the form 

AD qA (z) - AD qA (z) = N q z a (l - zf^ - z). (9) 

These parameters are also given in Table I. We also performed a fit using flavor symmetric 
fragmentation functions which we shall need for the discussion of Lambda production in 
pp collisions. The fit parameters are given in Table. II. In Fig. 2, we show the calculated 
fragmentation functions at Q 2 = M|. We note that our fragmentation functions describe the 
asymmetry in leading and non-leading particle production, as well as the lambda polarization 



measured in e + e annihilation at the Z pole, very well — as has been shown in Ref. [14] 



III. POLARIZED PROTON PROTON COLLISION 

Spin-dependent fragmentation of quarks can be studied in proton-proton collisions with 
one of the protons polarized ||. Here, many subprocesses may lead to the final state Lambda 
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so that one has to select certain kinematic regions to suppress the unwanted contributions. In 
particular, in order to test whether polarized up and down quarks do fragment into polarized 
Lambdas the rapidity of the produced Lambda has to be large, since at high rapidity, A's 
are mainly produced through valence up and down quarks. (We count positive rapidity in 
the direction of the polarized proton beam.) 

The difference of the cross sections to produce a Lambda with positive helicity through 
the scattering of a proton with positive/negative helicity on an unpolarized proton is given 
in leading order perturbative QCD (LO pQCD) by [] 

E C ^(AB ^C + X)= E C -^(A^B - + X) - E C ^(A^B - & + X) 
d A p c d s p c d s p c 

= J2j dx a dx b dz c Af Aa (x a , n^fsbixb, ^ADccizc, ij 2 ) 

abed 

s Ada . , r . ~ , . . 

— 2 — ^(ab -► cd)5(s + t + u). (10) 

Here, A/ J 4 a (x a , /i 2 ) and fBb{x b ,/J> 2 ) are the polarized and unpolarized distribution functions 
of partons a and b in protons A and B, respectively, at the scale fi. x a and x b are the 
corresponding momentum fractions carried by partons a and b. AD c c{z c , fi 2 ) is the polarized 
fragmentation function of parton c into baryon C, in our case C = A. z c is the momentum 
fraction of parton c carried by the produced Lambda. Ada/dt is the difference of the cross 
sections at the parton level between the two processes a) +b — > + d and er + b — > + d. 
The unpolarized cross section is given by Eq. ( 110) with the A's dropped throughout. 
The Mandelstam variables at the parton level are given by 

s = x a x b s, i = -x a p ± y/se~ y / z c , u = -x b p ± y r se y / z c (11) 

where, y and p± are the rapidity and transverse momentum of the produced Lambda and 
yfs is the total center of mass energy. The summation in Eq.(|TH|) runs over all possible 
parton-parton combinations, qq' — ► qq', qg — > qg, qq — > qq... The elementary unpolarized 
and polarized cross sections can be found in Refs. |24], 25| . Performing the integration in Eq. 
( PD over z c one obtains 



Ada r 1 r 1 

E C -^{AB -^C+X) = Y, dx a dx b Af Aa (x a , fi 2 )f Bb (x b , fi 2 )AD cC (z c , fi 2 ) 

1 Ada 



d 3 Pc 

abed Xarn ^ n x bmin 

(ab -> cd) (12) 



TTZ C dt 

with 



X± —y y X a X±C y ^ X±C y (TO\ 

2x~ b 6 + 2x~ a 6 ' Xbmin= 2x a -x ± ey 1 X ~ = 2 - x ± e~v [ ' 



Since the relevant spin dependent cross sections on the parton level are only known in LO we 
perform a LO calculation here. 
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where x± = 2p± / ^/s. 

In order to elucidate the kinematics, in Fig. §, we plotted z c as a function of x a and y 
for two different transverse momenta, p± = 10 GeV (left) and p± = 30 GeV (right) and 
for two different values of Xb, Xb = Xbmin + 0.01 (top) and and Xb = Xbmin + 0.1 (bottom) 
in Fig. [| Note, that z c is maximal both for Xb = Xbmin and x a = x amin . With increasing 
rapidity, y, both the lower integration limit of x a , x am i n , and the momentum fraction of 
the fragmenting quark transferred to the produced A, z c , increase. Hence, large rapidities 
probe the fragmentation of mostly valence quarks into fast Lambdas. The dependence 
on the transverse momenta is also shown in Fig. |^. With increasing p±, the kinematic 
boundary is shifted to smaller rapidities and the fragmentation of the valence quarks can 
be studied at lower rapidities. This is important since the available phase space is limited 
by the acceptance of the detectors at RHIC. However, the cross section also decreases with 
increasing transverse momenta, leading to lower statistics. 

In Fig. 4, we show the contributions of the various channels to the cross section for two 
different transverse momenta, both for inclusive Lambda (4a) and inclusive jet production 
(4b). gq —>■ gq stands, for example, for the contribution to the cross section coming from 
the subprocess involving a gluon g and a quark q in the initial and final states. In qql — > qq/, 
the quarks have different flavors and qq! — > qq! is also included. Although the kinematics 
are not exactly the same for these two processes one can study the role played by the 
fragmentation functions by comparing inclusive Lambda and inclusive jet production. In 
particular, the contributions from channels containing two gluons in the final state are 
suppressed in inclusive Lambda production due to the smallness of D g \. We note that 
D g \ has been set to zero at the initial scale and is generated through evolution. Thus, 
while qg — > qg is the dominant channel in inclusive Lambda production, both qg — > qg and 
99 ~~ * 99 are equally important in inclusive jet production. There is some ambiguity due to 
our poor knowledge of the gluon fragmentation — larger probabilities for g — > A will enhance 
the contributions from gluons in inclusive Lambda production. However, the contribution 
from the process gg — > gg falls off faster than that from qg — > qg with increasing rapidity, 
since g(x a ) decreases faster than q(x a ) with increasing x a and the Lambdas are produced 
mainly from valence up and down quarks at high rapidities. 

Our analysis of the kinematics and the various contributions to inclusive Lambda pro- 
duction already indicate that Lambda polarization measurements in pp collisions at high 
rapidities are ideally suited to test whether polarized up and down quarks may fragment 
into polarized A's. We calculated the Lambda polarization using our flavor asymmetric 
fragmentation functions for RHIC energies and for p± = 10 GeV. Increasing the transverse 
momentum gives similar results, with the only difference that Pa starts to increase at lower 
rapidities. We used the standard set of GRSV LO quark distributions for the polarized 
parton distributions |26| and the LO Cteq4 distributions for the unpolarized quark distribu- 
tions pTfl . The scale, fi, is set equal to p±. We also checked that there is only a very weak 
dependence on the scale by calculating the polarization using fi = p±/2 and fi = 2p±. The 



2 While there is only one integration variable, x a , in inclusive jet production, once p± and y are 
fixed, both x a and Xb have to be integrated over the allowed kinematic region in inclusive Lambda 
production, since the produced A carries only a fraction of the parton's momentum. 
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predicted Lambda polarization is shown in Fig |5]a. It is positive at large rapidities where 
the contributions of polarized up and down quarks dominates the production process. At 
smaller rapidities, where x a is small, strange quarks also contribute. However, since the 
ratios of the polarized to the unpolarized parton distributions are small at small x a the 
Lambda polarization is suppressed. The result also depends on the parameterization of the 
polarized quark distributions. In particular, the polarized gluon distribution is not well con- 
strained. However, it is clear from the kinematics that the ambiguity associated with the 
polarized gluon distributions only effects the results at lower rapidities. This can be seen 
in Fig. |5|b where we plot the contribution from gluons, up plus down quarks and strange 
quarks to the Lambda polarization. 

Next, we contrast our prediction with the predictions of various SU(3) flavor symmetric 
models which use 

D uA = D dA = D sA . (14) 

We fitted the cross sections in e + e~ annihilation using Eq. ([14]) and the functional form given 
in Eq. (^). For the polarized fragmentation functions, we discuss two different scenarios: 
The model, SU(3)a (c.f. Fig. 5a), corresponds to the expectations of the naive quark model 
that only polarized strange quarks can fragment into polarized Lambdas 

AD uA = AD dA = AD sA = D sA . (15) 

It gives essentially zero polarization because the strange quarks contribute at low rapidities 
where the polarization is suppressed. Model, SU(3)b (c.f. Fig. 5a), which was proposed in 
Ref. 0, is based on DIS data, and sets 

AD uA = AD dA = -0.20D uA AD sA = O.QOD sA . (16) 

This model predicts negative Lambda polarization. 

Finally, we address the problem of Lambdas produced through the decay of other hyper- 
ons, such as E° and E*. In order to estimate the contribution of hyperon decays we assume, 
in the following, that 

(1) the A's produced through hyperon decay inherit the momentum of the parent hyperon 

(2) and that the total probability to produce A, E° or E* from a certain uds state is 
given by the 577(6) wave function and is independent of the mass of the produced hyperon. 

Further, in order to estimate the polarization transfer in the decay process we use the 
constituent quark model. The polarization can be obtained by noting that the boson emitted 
in both the E° — > A7 and the E* — > An decay changes the angular momentum of the 
nonstrange diquark from J = 1 to J = 0, while the polarization of the spectator strange 
quark is unchanged. Then, the polarization of the A is determined by the polarization of the 
strange quark in the parent hyperon, since the polarization of the A is exclusively carried 
by the strange quark in the naive quark model. 

First, let us discuss the case when the parent hyperon is produced by a strange quark. 
Since the strange quark is always accompanied by a vector ud diquark, in both E° and E* 
the fragmenation functions of strange quarks into these hyperons are much softer than the 
corresponding fragmentation function into a A. Thus, in the high z limit, the contributions 
from the processes, s —>■ E° — > A and s —>■ E* — > A, are negligible compared to the direct 



S 



production, s — > A. Furthermore, both channels, s — > S° — > A and s — > X* — > A, enhance 
the already positive polarization from the direct channel, s — > A. 

This is different in the case when the parent hyperon is produced by an up or down quark. 
Both A and S° can be produced by an up (down) quark and a scalar ds (us) diquark — a 
process which dominates in the large z limit. (The component with a vector diquark can be 
neglected in this limit). Furthermore, the up and down fragmentation function of the S* are 
as important as those of the A and S° in the large z limit. This is because the u fragmentation 
function of £* peaks at about 1385/(1010 + 1385) ~ 0.58 which is almost the same as the 
peak of the scalar components of the A and S, which are 1115/(890 + 1115) ~ 0.57 and 
1190/(890 + 1190) ~ 0.57, respectively. Thus, for the up and down quark fragmentation, it 
is important to include the A's from these decay processes. 

The relevant probabilities to produce a A with positive and negative polarization from a 
fragmenting up quark with positive polarization and an ds diquark are shown in Table |T|. 
We assumed that all spin states of the ds diquark are produced with equal probabilities. 
The final weights which are relevant in the large z limit are set in bold. We find that if we 
include all channels, which survive in the large z limit, the polarization of the A is reduced 
by a factor of 10/27 compared to the case where only the directly produced A's are included. 
Since the £* decay is a strong decay it is sometimes included in the fragmentation function 
of the A. Including only £*, the suppression factor we obtain is 49/81. (Note that our model 
predicts that u^(ds) 0i o — > A, u^(ds) fl — ► S° and u^(ds)i >0 — » £* have approximately the 
same z dependence and are approximately equal (up to the Clebsch- Gordon factors) since 
the ratios, M/(M + M n ), have roughly the same numerical values. Thus, the effect of the 
E° and £* decays can be taken into account by a multiplicative factor.) 

In order to illustrate the effect of these decays on the final A polarization, we multiplied 
our results with these factors. The results are shown in Fig. 6b as dotted lines. We note 
that our implementation of this correction relies on the assumptions that the produced A 
carries all the momentum of the parent hyperon and that all states are produced with equal 
probabilities. Since neither of these assumptions is strictly valid, we tend to overestimate 
the importance of hyperon decays. Note also that the inclusion of S° decay in the £77(3) 
symmetric models makes the resulting polarization more negative. As a result, even if effects 
of £ decays are included, large discrepancies still persist between our predictions and those 
of ££7(3) symmetric models. 

IV. CONCLUSIONS 

Measurements of the Lambda polarization at RHIC would provide a clear answer to the 
question of whether polarized up and down quarks can transfer polarization to the final state 
A. We predict positive Lambda polarization at high rapidities, in contrast with models based 
on ££7(3) flavor symmetry and DIS which predict zero or negative Lambda polarization. Our 
prediction is based on the same physics which led to harder up than down quark distributions 
in the proton and to the A-N and S-A mass splittings. We also estimated the importance 
of S° and £* decays which tend to reduce the predicted A polarization. 
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TABLES 



TABLE I. Fit parameters obtained by fitting the e + e data. We also parametrized D q A — DqA 
and AD q A — ADqA, calculated in the bag. 



Parameter 


D S A - DgA 


D U A - D U A 


Dq A 


A -D s a - AD S - A 


AD uA - AD U A 


N 


5.81 x 10 9 


1.60 x 10 17 


99.76 


3.73 x 10 18 


-6.25 x 10 10 


a 


21.55 


30.49 


1.25 


21.21 


32.48 


(3 


13.60 


28.34 


11.60 


13.38 


27.72 


7 










0.52 



TABLE II. Fit parameters obtained by fitting the e + e data and asumming that the fragmen- 
tation fucntions are flavor symmetric. 
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DqA - DqA 


DqA 
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1.92 x 10 4 


99.76 
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7.47 


1.25 
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8.06 


11.60 
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TABLE III. Different channels for the production of A hyperons from a positively polarized up 
quark and a ds diquark. It is assumed that all spin states of the ds diquark are produced with the 
same probabilities. S* T and stand for the 1/2 and 3/2 spin component of the S*. See text for 
further details. 
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FIG. 1. Inclusive Lambda production in e + e~ annihilation. The solid lines are the result of the 
global fit. They contain two parts, the fixed contributions from D q \ — Dg\ calculated in the bag 
(dashed line only shown for the Aleph data) and Dg\ obtained from the fit (dash-dotted line), xe 
is defined as xe = 2E\/y/s where E\ is the energy of the produced Lambda in the e + e~ center of 
mass frame and ^/s is the total center of mass enery. 
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FIG. 2. Fragmentation functions. The solid and dashed lines stand for the calculated fragmen- 
tation functions of up and strange quarks into Lambda baryons through production of a Lambda 
and an anti-diquark and correspond to D u \ — D^a and D s \ — D s \, respectively. The dash-dotted 
line represents the contributions from higher intermediate states, and is obtained by fitting the 
e + e~ data and corresponds to Dq\. The short dashed line is the gluon fragmentation function. 
The light and heavy lines are the fragmentation functions at the scales Q 2 = /i 2 and Q 2 = Mf, 
respectively. Note that D g \ = at Q 2 = fi 2 . 
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FIG. 4. Contributions from the various channels (a) to the inclusive Lambda production cross 
section {pp — ► A + X) and (b) to the inclusive jet production cross section (pp — ► jet + X) at 
p ± = 10 GeV (left) and p ± = 30 GeV (right) at ^ = 500 GeV. 
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FIG. 5. Lambda polarization at RHIC. (a) The solid line represents our prediction. The pre- 
dictions of SU(3) symmetric fragmentation models are shown for comparision. The model labeled 
as SU(3)a is based on the quark model expectation that only the polarized strange quark may 
fragment into polarized Lambdas, while SU(3)b, is based on DIS data, (b) Contributions of dif- 
ferent flavors to the A-polarization. The light dashed, dash-dotted and heavy dashed lines stand 
for the contributions from up plus down, from strange and from gluon fragmentation, respectively, 
as calculated here. The estimated polarization including both S° and E* (lower dotted line) and 
only E* (upper dotted line) decays are also shown. See text for further details. 
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